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Abstract 

We present some phenomenology of a new class of intersecting D-brane models. Soft 
SUSY breaking terms for these models are calculated in the it-moduli dominant SUSY 
breaking approach (in type 1 1 A). In this case, the dependence of the soft terms on the 
Yukawas and Wilson lines drops out. These soft terms have a different pattern compared 
to the usual heterotic string models. Phenomenological implications for dark matter are 
discussed. 
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1 Introduction 



One of the goals of string phenomenology is to explain/predict features of low energy physics - 
both qualitatively and quantitatively. We are still far from that elusive goal. To make progress 
we think it is essential to build more and more realistic string models and to study their 
phenomenological features. Until a few years ago, only the heterotic string was considered a 
serious candidate for providing the unified theory of fundamental interactions. For a sample 
of heterotic string model building, see pQ. Developments in the past few years have shown 
that type I and type II strings provide us with new classes of J\f — 1, D — 4 vacua, with 
new avenues for model building. In addition, the concept of D-branes has provided us with 
a better understanding of type / (or equivalently type IIB orientifold) string theory. It has 
recently become evident that intersecting D-brane models offer excellent opportunities for string 
phenomenology. In fact, these developments have been collectively dubbed as the second string 
(phenomenology) revolution 

This paper is devoted to the detailed study of a particular class of models based on type II 
string theory compactifications on Calabi-Yau manifolds with Dj?-branes wrapping intersecting 
cycles on the compact space. This approach to string model building is distinguished by its 
computability and simplicity, together with very appealing phenomenological possibilities. In 
these models, gauge interactions are confined to D-branes. Chiral fermions are open strings 
which are stretched between two intersecting branes. They are localized at the brane inter- 
sections. If certain conditions are satisfied, there will be massless scalars associated with the 
chiral fermions such that we have Af = 1 supersymmetry in the effective field theory. Because 
of these attractive features, intersecting brane model building has drawn considerable attention 
in recent years and several semi-realistic models with an SM or MSSM like spectrum have been 
constructed |3J El ■ 

To test these approximate clues and to begin to test string theory, only reproducing the SM 
particle content is not enough. Numerical predictions must be made. In addition, a successful 
theory should not just explain existing data, it must also make predictions which can be tested 
in future experiments. For the brane models, if supersymmetry exists and is softly broken, 
soft SUSY breaking terms can calculated and tested by future experimental measurements. 
A fair amount of work on the low-energy effective action of intersecting D-brane models has 
been done. The stability of these kind of models has been discussed in The question of 
tree level gauge couplings, gauge threshold corrections and gauge coupling unification has been 
addressed in 13 UH HI] ■ Yukawa couplings and higher point scattering have been studied 
in ^21 E21 El E3- Some preliminary results for the Kahler metric have been obtained in [TE| . 
A more complete derivation of the Kahler metric directly from open/closed string scattering 
amplitudes has been done in (17] . which we use in this paper. The question of supersymmetry 
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breaking has also been addressed in such models [TBI dj 1201 EH 122] , using techniques of flux 
compactification in type I IB string theory. The basic idea is that turning on background closed 
string 3-form fluxes can provide a source of supersymmetry breaking. Although an exciting 
idea, it is not yet complete in all respects, see [23] . 

In this paper, we have taken a more phenomenological approach, parametrizing the effects of 
supersymmetry breaking in a more model independent manner and examining the consequences. 
Our main goal here is to use the results of [T7j to calculate and analyze effective low energy 
soft supersymmetry breaking terms. We also look at some of their dark matter applications. 
Applications to collider phenomenology will be dealt with in future work. Our main purpose in 
this paper is to move the string constructions of this approach closer to broad phenomenological 
applications. While we were writing the manuscript, [20] appeared, which has a large overlap 
in the computation of the soft terms. 

The paper is organized as follows. In section |2J we briefly review the intersecting brane 
model constructions. Then in section El we describe a brane setup in detail, for which we 
will study the soft terms. This brane setup was first introduced in [T2] . We compute the 
soft terms under the assumption of section HJ in the w-moduli SUSY breaking scenario. We 
compute the general formulas for the soft terms so that they can be applied to a large class of 
type IIA/B brane setups, including the flux compactification approach [TBI HOI 1201 I2JJ 122] • To 
understand opportunities better, we apply these general formulas to a particular brane setup, 
and study three particular points in the parameter space. One of these gives a W LSP, the 
second gives a H LSP and the third gives a mixed B-H LSP. The three points represent an 
almost generic feature of the parameter space of these intersecting brane models if one requires 
a light gluino, which amounts to reducing fine-tuning as suggested in [24J. In section 03 we 
discuss the phenomenological implications of the above model: the structure of soft terms, 
spectrum, gauge unification, issues of flavor and phase and in particular, the consequences for 
cosmology. We conclude in section |01 Some technical details are provided in the Appendix. 

2 General construction of intersecting brane models. 

In this section, we will briefly review the basics of constructing these models. More comprehen- 
sive treatments can be found in [2H1 1201 12H I2H1 120] ■ The setup is as follows - we consider type 
II A string theory compactified on a six dimensional manifold M.. It is understood that we are 
looking at the large volume limit of compactification, so that perturbation theory is valid. In 
general, there are K stacks of intersecting D6-branes filling four dimensional Minkowski space- 
time and wrapping internal homology 3-cycles of A4. Each stack P consists of Np coincident 
D6 branes whose worldvolume is M 4 x Up, where Up is the corresponding homology class of 
each 3-cycle. The closed string degrees of freedom reside in the entire ten dimensional space, 
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which contain the geometric scalar moduli fields of the internal space besides the gravitational 
fields. The open string degrees of freedom give rise to the gauge theory on the D6-brane world- 
volumes, with gauge group lip U(Np). In addition, there are open string modes which split into 
states with both ends on the same stack of branes as well as those connecting different stacks 
of branes. The latter are particularly interesting. If for example, the 3-cycles of two different 
stacks, say lip and Uq intersect at a single point in A4, the lowest open string mode in the 
Ramond sector corresponds to a chiral fermion localized at the four dimensional intersection of 
P and Q transforming in the bifundamental of U(Np) x U(Nq) [3U]. The net number of left 
handed chiral fermions in the ab sector is given by the intersection number Ipq = [Hp] ■ [IIq]. 

The propagation of massless closed string RR modes on the compact space M. under which 
the D-branes are charged, requires some consistency conditions to be fulfilled. These are known 
as the RR tadpole-cancellation conditions, which basically means that the net RR charge of 
the configuration has to vanish [31]. In general, there could be additional RR sources such as 
orientifold planes or background fluxes. So they have to be taken into account too. Another 
desirable constraint which the models should satisfy is M = 1 supersymmetry. Imposing this 
constraint on the closed string sector requires that the internal manifold M. be a Calabi-Yau 
manifold. We will see shortly that imposing the same constraint on the open string sector leads 
to a different condition. 

A technical remark on the practical formulation of these models is in order. Till now, we 
have described the construction in type II A string theory. However, it is also possible to 
rephrase the construction in terms of type II B string theory. The two pictures are related 
by T-duality. The more intuitive picture of type II A intersecting D-branes is converted to a 
picture with type IIB D-branes having background magnetic fluxes on their world volume. 
It is useful to remember this equivalence as it turns out that in many situations, it is more 
convenient to do calculations in type IIB. 

Most of the realistic models constructed in the literature involve toroidal (T 6 ) compacti- 
fications or orbifold/orientifold quotients of those. In particular, orientifolding introduces 06 
planes as well as mirror branes wrapping 3-cycles which are related to those of the original 
branes by the orientifold action. For simplicity, the torus (T 6 ) is assumed to be factorized into 
three 2-tori, i.e T 6 = T 2 x T 2 x T 2 . Many examples of the above type are known, especially 
with those involving orbifold groups - i) Z 2 x Z 2 |H2] ii) Z± x Z 2 [33], hi) Z4 EH], iv) Z G [33] . 
etc. 



3 A local MSSM-like model 

In order to make contact with realistic low energy physics while keeping supersymmetry intact, 
we are led to consider models which give rise to the chiral spectrum of the MSSM. It has been 
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Stack 


Number of Branes 


Gauge Group 




(nl, ml) 


(n 3 a ,m 3 a ) 


Baryonic 


N a = 3 


U(3) = SU{3) x U(l) a 


(1,0) 


(1/P,3p) 


(1/P, -3p) 


Left 


N b = l 


USp{2) = SU{2) 


(0,1) 


(1,0) 


(o,-i) 


Right 


N c = l 


U(l) e 


(0,1) 


(o,-i) 


(1,0) 


Leptonic 


N d =l 


U(l) d 


(1,0) 


(1/P,3p) 


(VP, -3p) 



Table 1: Brane content for an MSSM-like spectrum. The mirror branes a*,b*,c*,d* are not shown. 
p can take values 1, f/3. For concreteness, we take p = 1 for calculating the soft terms. However, the 
parameter space for the soft terms remains the same for both p = 1 and p = 1/3. 

shown in |3] that this requires us to perform an orientifold twist. A stack of Np D6 branes 
wrapping a 3-cycle not invariant under the orientifold projection will yield a U(Np) gauge 
group, otherwise we get a real {SO{2Np)) or pseudoreal {USp {2Np)) gauge group. 

Using the above fact, the brane content for an MSSM-like chiral spectrum with the correct 
intersection numbers has been presented in |12| . Constructions with more than four stacks of 
branes can be found in [HE]- In the simplest case, there are four stacks of branes which give 
rise to the initial gauge group : U{3) a x Sp{2)b x U{l) c x U{l)d, where a,b,c&:d label the 
different stacks. The intersection numbers Ipg = [Hp] ■ [IIq] between a D6-brane stack P and 
a D6-brane stack Q is given in terms of the 3-cycles [lip] and [IIq], which are assumed to be 
factorizable. 

[lip] = [{np, nip) ® (n P , mp) <g> {np, mp)} (1) 

where (n P , m P ) denote the wrapping numbers on the i th 2-torus.The 06 planes are wrapped 
on 3-cycles : 

[n O6 ] = (g)[(i,0)r (2) 

r=l 

Note that for stack b, the mirror brane b* lies on top of b. So even though Nf, = 1, it can 
be thought of as a stack of two D6 branes, which give an USp{2) = SU{2) group under the 
orientifold projection. 

The brane wrapping numbers are shown in Table ^ and the chiral particle spectrum from 
these intersecting branes are shown in Table 121 

3.1 Getting the MSSM 

The above spectrum is free of chiral anomalies. However, it has an anomalous U{1) given by 
U{l) a + U{l)d- This anomaly is canceled by a generalized Green-Schwarz mechanism [26], 
which gives a Stuckelberg mass to the U{1) gauge boson. The two nonanomalous U{l)s are 
identified with {B — L) and the third component of right-handed weak isospin U{l)p [12] . In 
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fields 


sector 


I 


SU(3) C x SU(2) L 


U(l)a 


U{l)c 


u(i) d 


U(l)y 


Ql 


(a,b) 


3 


(3,2) 


1 








1/6 


Ur 


(a,c) 


3 


(3,1) 


-1 


1 





-2/3 


Dr 


(a,c*) 


3 


(3,1) 


-1 


-1 





1/3 


L 


(d,b) 


3 


(1,2) 








1 


-1/2 


E R 


(d,c*) 


3 


(1,1) 





-1 


-1 


1 


Nr 


(d,c) 


3 


(1,1) 





1 


-1 





H u 


(b,c) 


1 


(1,2) 





-1 





1/2 


H d 


(b,c*) 


1 


(1,2) 





1 





-1/2 



Table 2: The MSSM spectrum from intersecting branes. The hypercharge normalization is 
given by Q Y = \Q a - \Qc - \Qd- 

orientifold models, it could sometimes happen that some nonanomalous U(l)s also get a mass 
by the same mechanism the details of which depend on the specific wrapping numbers. It 
turns out that in the above model, two massless U(l)s survive. In order to break the two U{l)s 
down to U(1)y, some fields carrying non- vanishing lepton number but neutral under U(l)y 
are assumed to develop vevs. This can also be thought of as the geometrical process of brane 
recombination [2TJ EH] ■ 

3.2 Global embedding and supersymmetry breaking 

As can be checked from Table 1, the brane content by itself does not satisfy the RR tadpole 
cancellation conditions : 

£([n a ] + [n Q .]) =32[n 06 ] (3) 

a 

Therefore, this construction has to be embedded in a bigger one, with extra RR sources in- 
cluded. There are various ways to do this such as including hidden D-branes or adding back- 
ground closed string fluxes in addition to the open string ones. As a bonus, this could also 
give rise to spontaneous supersymmetry breaking. With extra D-branes, one might consider 
the possibility of gaugino condensation in the hidden sector [3HJ. Alternatively, one could con- 
sider turning on background closed string NS-NS and RR fluxes which generate a non-trivial 
effective superpotential for moduli, thereby stabilizing many of them [T%1 I19| 12*2*]. 

In this paper, we will leave open the questions of actually embedding the above model 
in a global one and the mechanism of supersymmetry breaking. We shall assume that the 
embedding has been done and also only parametrize the supersymmetry breaking, in the spirit 
of |3*9*1 HP] . We are encouraged because there exists a claim of a concrete mechanism for the 
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global embedding of (the T-dual of) this model as well as supersymmetry breaking [2T] . 
3.3 Exotic matter and fi problem 

The above local model is very simple in many respects, especially with regard to gauge groups 
and chiral matter. However, it also contains exotic matter content which is non-chiral. These 
non-chiral fields are related to the untwisted open string moduli - the D-brane positions and 
Wilson lines. The presence of these non-chiral fields is just another manifestation of the old 
moduli problem of supersymmetric string vacua. However, it has been argued [201 IS] that 
mass terms for the above moduli can be generated by turning on a F- theory 4-form flux. One 
then expects that a proper understanding of this problem will result in a stabilization of all 
the moduli. As explained in |21j . there could be Af = 1 embeddings of this local model in 
a global construction. This requires additional D-brane sectors and background closed string 
3-form fluxes. The other D-brane sectors add new gauge groups as well as chiral matter, some 
of which could be charged under the MSSM gauge group. This may introduce chiral exotics in 
the spectrum, an undesirable situation. However, many of these exotics uncharged under the 
MSSM gauge group can be made to go away by giving vevs to scalars parametrizing different 
flat directions. In this paper, we assume that there exists an embedding such that there are no 
chiral exotics charged under the MSSM. Such exotics can cause two types of problems. It is 
of course essential that no states exist that would already have been observed. It seems likely 
that can be arranged. In addition, states that would change the RGE running and details of 
the calculations have to be taken into account eventually. 

The higgs sector in the local model arises from strings stretching between stacks b and c. 
However, the net chirality of the be sector is zero, since the intersection number l\, c is zero. The 
higgs sector in the above model has a \i term, which has a geometrical interpretation. The 
real part of the \i parameter corresponds to the separation between stacks b and c in the first 
torus, while the imaginary part corresponds to a Wilson line phase along the 1-cycle wrapped 
on the first torus. These correspond to flat directions of the moduli space. Adding background 
closed string fluxes may provide another source of /i term JH] > which will lift the flat direction 
in general. Thus, the effective \i term relevant for phenomenology is determined by the above 
factors and the problem of obtaining an electroweak scale /i term from a fundamental model 
remains open. In this paper, therefore, we will not attempt to calculate /z, and fix it by imposing 
electroweak symmetry breaking (EWSB). It is important to study further the combined effect 
of the several contributions to /i and to EWSB. 
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3.4 Type II A - type IIB equivalence 

As mentioned earlier, it is useful to think about this model in terms of its T-dual equiva- 
lent. In type IIB, we are dealing with D9 branes wrapped on T 2 x T 2 x T 2 with an open 
string background magnetic flux T 3 turned on. Therefore the D9-branes have in general mixed 
Dirichlet and Neumann boundary conditions. The flux has two parts - one coming from the 
antisymmetric tensor (b 3 ) and the other from the gauge flux (F J ) so that : 

P = V + 2na' F j (4) 

The above compactification leads to the following closed string Kahler and complex structure 
moduli, each of which are three in number for this model: 



R 



T'i = V + %R(R 3 sin(cO; U' j = Hj e ia " ; j = l, 2, 3. (5) 

R 1 

where i?/ and R^ are lengths of the basis lattice vectors characterizing the torus T 2,3 and a? 
is the angle between the two basis vectors of the torus T 2j . By performing a T-duality in the 
y direction of each torus T 2j , the D9 brane with flux T 3 is converted to a D6 brane with an 
angle Qi with respect to the x-axis. This is given by [42J: 

where fi is defined by the quantization condition for the net 2-form fluxes P as 

Using the above equation and the relation between the type II A and type IIB: 

a' ■ a' 

T 3 = -tJt, U 3 = -^r, J = 1,2,3. (8) 



R{ - Rf, R\ - —77 (9) 

R-2 



we get the corresponding type II A relation: 

m? \U ] \ 2 

^> = wikm ,i = w (10) 

The unprimed symbols correspond to the type IIA version while the primed ones to the type 
IIB. 
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3.5 Af= 1 SUSY 



We now look at the M = 1 supersymmetry constraint on the open string sector. In type 1 1 A, 
this leads to a condition on the angles 9^ [30 : 

6 1 + 6 2 + 6 3 = mod 2 or 

3 m J I TP 1 2 3 I TP 1 2 

V— ]U 1 = TT — ' ' (11) 
7^' Im (Etf ) ni Im (t/J) 1 ; 

which after T-duality leads to a condition on the fluxes in type IIB. 



4 Low energy effective action and soft terms 



We now analyze the issue of deriving the four dimensional M = 1 low energy effective action 
of these intersecting brane models. In the type IIB picture, this has been done in |43t l4*4"j by 
Kaluza Klein reduction of the Dirac-Born-Infeld and Chern-Simons action. The effective action 
can also be obtained by explicitly computing disk scattering amplitudes involving open string 
gauge and matter fields as well as closed string moduli fields and deducing the relevant parts 
of the effective action directly. This has been done in ^7]. We will follow the results of ^7j in 
our analysis. 

The J\f = 1 supergravity action thus obtained is encoded by three functions, the Kahler 
potential K, the superpotential W and the gauge kinetic function / ■ Each of them will 
depend on the moduli fields describing the background of the model. One point needs to be 
emphasized. When we construct the effective action and its dependence on the moduli fields, 
we need to do so in terms of the moduli s, t and u in the field theory basis, in contrast to the 
S, T and U moduli in the string theory basis J7j. In type 1 1 A, the real part of the field theory 
s, u and t moduli are related to the corresponding string theory S, U and T moduli by : 



Re(s) 



Vim U 1 Im U 2 Im U r 



2tt 



2tt 



Im U k Im U l 



(j,k,l) = (1,2,3) 



lot 



(12) 



where j stands for the j th 2-torus. The above formulas can be inverted to yield the string 
theory U moduli in terms of the field theory moduli s and u. 



\U j \ 



Re(w fc )Re(^ 



lm(W) \ Re(vP)Re(s) 



(j,k,l) = (1,2,3) 



(13) 
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The holomorphic gauge kinetic function for a D6 brane stack P is given by : 

e 1 / 1 2 S 12 3 1 2 1^2 S 1 2 S\ 

j P = — (rip n p n p s — n p m p m p u — n p m p m p u — n P m p m p u ) 

Kp 

9m P = |Re(/p)| (14) 

The extra factor is related to the difference between the gauge couplings for U(Np) and 
Sp(2N P ), SO(2N P ). k p = 1 for U(N P ) and k p = 2 for Sp(2N P ) or SO(2N P ) 
The SM hypercharge U(l)y gauge group is a linear combination of several U(l)s: 

Qy = l -Q a - \q c - l -Q d - (15) 

O 2 I 

Therefore the gauge kinetic function for the U(l)y gauge group is determined to be[TT]: 

fy = g/D a + -Id c + -^fD d - (16) 
The Kahler potential to the second order in open string matter fields is given by : 
K(M,M,C,C)=K(M,M)+ K CiC] {M,M)C % C 3 + £ K CgCe (M,M)C e C e (17) 

untwisted i,j twisted, 9 

where M collectively denote the moduli; Cj denote untwisted open string moduli which comprise 
the D-brane positions and the Wilson line moduli which arise from strings with both ends on 
the same stack; and Cg denote twisted open string states arising from strings stretching between 
different stacks. 

The open string moduli fields could be thought of as matter fields from the low energy field 
theory point of view. The untwisted open string moduli represent non-chiral matter fields and 
so do not correspond to the MSSM. For the model to be realistic, they have to acquire large 
masses by some additional mechanism, as already explained in section 

Let's now write the Kahler metric for the twisted moduli arising from strings stretching 
between stacks P and Q, and comprising 1/4 BPS brane configurations. In the type IIA 
picture, this is given by 12313 UZj : 

where 6 P q — Op — On is the angle between branes in the j th torus and 04 is the four dimensional 
dilaton. From (fEZj) . 04 can be written as (Re(s) Re(iti) Re(ii 2 ) Re(w 3 ))~ 1//4 . The above Kahler 
metric depends on the field theory dilaton and complex structure moduli u 1 through 4 and 
9pq. It is to be noted that (j!8)l is a product of two factors, one which explicitly depends on 



rn - o\ 



PQ) 



r(oi 



PQ) 
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the field theory s and u moduli (e^ 4 ), and the other which implicitly depends on the s and u 
moduli (through the dependence on 6 P q). Thus, Kpg can be symbolically written as : 

K PQ = e* K° PQ (19) 

The Kahler metric for 1/2 BPS brane configurations is qualitatively different from that of the 
1/4 BPS brane configurations mentioned above. Generically, these cases arise if both branes 
P and Q have a relative angle 8 P q = 0, 1 in the same complex plane j. It is worthwhile to 
note that the higgs fields in Table El form a 1/2 BPS configuration and are characterized by the 
following Kahler metric [T9*] : 

K^ s =((s + s)(u 1 + uW + m 3 + ?)y 1/2 (20) 

An important point about the above expressions needs to be kept in mind. These expressions 
for the Kahler metric has been derived using tree level scattering amplitudes and with Wilson 
line moduli turned off. Carefully taking the Wilson lines into account as in |13J, we see that the 
Kahler metric has another multiplicative factor which depends on the Wilson line moduli as 
well as t moduli in type II A. If the Wilson line moduli do not get a vev, then our analysis goes 
through. However, if they do, they change the dependence of the metric on the t moduli. As 
will be explained later, we only choose the u moduli dominated case for our phenomenological 
analysis, so none of our results will be modified. 
The superpotential is given by: 

W = W + ~/i a/3 (M) C a + i Y af3l {M) C a C 3 C + ... (21) 

In our phenomenological analysis, we have not included the Yukawa couplings for simplicity. 
But as explained later, in the u moduli dominant SUSY breaking case, the soft terms are 
independent of the Yukawa couplings and will not change the phenomenology. 

4.1 Soft terms in general soft broken M = 1, D = 4 supergravity 
Lagrangian 

From the gauge kinetic function, Kahler potential and the superpotential, it is possible to 
find formulas for the normalized soft parameters - the gaugino mass parameters, mass squared 
parameter for scalars and the trilinear parameters respectively. These are given by [4*U] : 




m 2 PQ = (mj /2 +Vo)-£ F^F N d Al d N log(K PQ ) 

M,N 

A PQR = F M [k M + d M \og(Y PQR )-d M \og(kp Q K QR K RP )] (22) 
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For our purposes, P, Q and R denote brane stacks. So Mp denotes the gaugino mass parameter 
arising from stack P; m P Q denotes mass squared parameters arising from strings stretching 
between stacks P and Q and Ypqp denotes Yukawa terms arising from the triple intersection 
of stacks P, Q and R. The terms on the RHS without the indices P, Q and R are flavor 
independent. Also, M and N run over the closed string moduli. F M stands for the auxiliary 
fields of the moduli in general. Supersymmetry is spontaneously broken if these fields get non- 
vanishing vevs. It is assumed here that the auxiliary fields D have vanishing vevs. Their effect 
on the soft terms can be calculated as in $T\, which we assume to be zero. These formulas 
have been derived for the case when the Kahler metric for the observable (MSSM) fields is 
diagonal and universal in flavor space. In principle, there are also off-diagonal terms in the 
Kahler metric. They relate twisted open string states at different intersections and therefore are 
highly suppressed. We neglect the off- diagonal terms in our study. If the seperations between 
the intersections are very small, the off-diagonal terms or non-universal diagonal terms may 
have observable effects leading to interesting flavor physics. 

The effective M = 1, 4d field theory is assumed to be valid at some high scale, presumably 
the string scale. The string scale in our analysis is taken to be the unification scale (~ 2 x 
10 16 GeV). We then need to use the renormalization group equations (RGE) to evaluate these 
parameters at the electroweak scale. In this paper, as mentioned before, it is assumed that the 
non-chiral exotics have been made heavy by some mechanism and there are no extra matter 
fields at any scale between the electroweak scale and the unification scale. This is also suggested 
by gauge coupling unification at the unification scale. 

One might wonder whether including the Yukawas in the analysis may lead to significant 
modifications in the spectrum at low energies because of their presence in the formulas for the 
soft terms (|22j). However, this does not happen. This is because the Yukawa couplings (Y a/ 3 7 ) 
appearing in the soft terms are not the physical Yukawa couplings (Y^S). The two are related 
by: 

W* 

Y$r = Y aPl ^ e*/ 2 (K a KpK 7 y^ (23) 

The Yukawa couplings (Y aj g 7 ) between fields living at brane intersections in intersecting D-brane 
models arise from worldsheet instantons involving three different boundary conditions j2Sj- 
These semi-classical instanton amplitudes are proportional to e~ A where A is the worldsheet 
area. They have been shown to depend on the Kahler (t) moduli (complexified area) and the 
Wilson line moduli [T2| in type II A. Although the physical Yukawas (Y^g ) depend on the u 
moduli through their dependence on the Kahler potential, the fact that Y" Q/ g 7 do not depend on 
the u moduli in type II A ensures that in the F u dominant supersymmetry breaking case, the 
soft terms are independent of Y a/3 ^. 

Thus our analysis is similar in spirit to those in the case of the heterotic string, where 
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dilaton dominated supersymmetry breaking and moduli dominated supersymmetry breaking 
are analyzed as extreme cases. It should be remembered however, that T-duality interchanges 
the field theory (t) and (u) moduli. Thus what we call u moduli in type II A, become t moduli 
in type I IB and vice versa. In a general situation, in which the F-terms of all the moduli get 
vevs, the situation is much more complicated and a more general analysis needs to be done. 
This is left for the future. 



4.2 Soft terms in intersecting brane models 

We calculate the soft terms in type IIA picture. As already explained, we assume that only 
F-terms for u moduli get non-vanishing vevs. In terms of the goldstino angle 9 as defined in 
|%0j . we have 9 = 0. We assume the cosmological constant is zero and introduce the following 
parameterization of F u *: 

F ui = v / 3m 3/2 (?i' + for i = 1, 2, 3 (24) 

with J2 |@«| 2 = 1- To calculate the soft terms, we need to know the derivatives of the Kahler 
potential with respect to u. For a chiral field C arising from open strings stretched between 
stacks P and Q, we denote its Kahler potential as Kpq. From fi"9|) . we see that their derivatives 
with respect to u % are 

d log Kpq = * d log K PQ d9 PQ -1 

du { py depQ du 1 4(V + m 1 ) 1 ' 

d 2 \ogk PQ = 3/ a log Kpg d 2 9 k PQ d 2 log Kpq dd k PQ dd k PQ 8jj \ 

dvtdui ^ I de k PQ dvtdv? d{9 k PQ ) 2 du { dti A^ + u 1 ) 2 ) [ ' 

From the Kahler potential in eq. (|18|) . we have 

dO J PQ 2 d6 PQ 2 d6 PQ 



d{o PQ y ' de j PQ 

The angle 9 P q = 9 p — 9q can be written in terms of u moduli as: 



[PQ) d(9 PQ y d9 PQ 1 



rrii 



I Reu k Reu l 



^^ffe* whereat,/) = (1,2,3) (29) 
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Then we have 



r)f) 3 

{u k + u k ) ^ 
' du k 



1 f 

j- sin(27r# J ') when j = k 



1 * 

sin(27T0 J ) when j ^ k 



(30) 



where we have defined [f(0^)]n = f{d 3 p) — f{&h)- The second order derivatives are 



0£g = (?/+u fc )(u ( +u 



T ^[sin(4vr^)+4sin(2vr^)]J 
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^[sin(47r^)-4sin(27r^)]J 



du k du l 



[sin(47r^)]J 



when j = k = I 

when j ^ k = I 

when j = ky^loTj = l^k 
when j 7^ k I j 



(31) 

Now, one can substitute eq. (|24H3lj) to the general formulas of soft terms in eq.(|22|). The 
formulas for the soft parameters in terms of the wrapping numbers of a general intersecting 
D-brane model are given by: 



Gaugino mass parameters 



Re/p i=1 v 



ripmpmp 



(j,k,l) = (1,2,3) 



(32) 



As an example, the can be obtained by putting P = a and using the appropriate 
wrapping numbers, as in Table H 



Trilinears : 

3 

Apqr = -v / 3m 3 / 2 51 



i + (E^(g-J) + (E^L. 4 , 



40 



fe=i 



fc=i 



+ 



2 



— m 3/2 Qie 471 



(33) 



This arises in general from the triple intersections PQ, QR and RP, where PQ and RP 
are 1/4 BPS sector states and QR is a 1/2 BPS state. QR being the higgs field, has a 
special contribution to the trilinear term compared to PQ and RP. So as an example, 
Aq l h u u r can be obtained as a triple intersection of ab, be and ca, as seen from Table El 
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Scalar mass squared (1/4 BPS) 



2 2 
m PQ = m 3/2 



m,n=l V j=l 



(3 

By using the definitions in eq. (j3T?l ETJ), we see that 9pQ and 6*pq do not depend on the 
vevs of u moduli: v? . As an example, the squark mass squared mj~ can be obtained by 
putting P,Q = a, b; as can be seen again from Tabled 

We can now use the wrapping numbers in Tabled to get explicit expressions for the soft terms 
for the particular model in terms of the s, t and u moduli and the parameters 7713/2, 7i and Oj, 
i = 1, 2, 3. The expressions for the trilinear A parameters and scalar mass squared parameters 
(except those for the up and down type higgs) are provided in the Appendix. Using Table HJ 
the formula for the gaugino mass parameters and the mass squared parameters for the up and 
down higgs are given by: 

• Gaugino mass parameters: 

3/2 3Rew 3 + 4Res + 36Re Ml v ; 

= V3m 3/2 e-^0 2 (36) 
R 9e~^ Rem Q 1 

M ° = m3//2 Reg + 9Rewi (37) 

It is important to note that there is no gaugino mass degeneracy at the unification scale, 
unlike other models such as mSUGRA. This will lead to interesting experimental signa- 
tures. 



Higgs mass parameters 



rnl u =rn% d = ml /2 {l-^\Q x \ 2 ) (38) 



The Higgs mass parameters are equal because they are characterized by the same Kahler 
metric Kp^ gs in (UJ). 

For completeness, we would also like to list the relative angles between a brane P and the 
orientifold plane on the jth torus are denoted by 9 P in Table El The soft terms depend on these 
angles. 
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Stack 


Number of Branes 


Gauge Group 


(o P ) 


(0%) 


(Op) 


Baryonic 


N a = 3 


U{3)=SU{3) x U(l) a 





a 


—a 


Left 


N b = l 


USp{2) = SU(2) 


1 

2 





l 

2 


Right 


N c = l 


U(l) c 


1 

2 


i 

2 





Leptonic 


N d = l 







a 


—a 



Table 3: Angles between a brane P and the orientifold plane on the jth torus: 9 3 P . 
The only non-trivial angle is a, which is given by 

^™>= s Wili- ( 3 9 > 

as can be seen from eq. f)10!13|) and Table [TJ 

We would now like to compare our setup and results for the soft terms with those obtained 
in [20] • The setup considered in j2Q| is very similar to ours, but with a few differences. It 
is a type I IB setup with three stacks of D7 branes wrapping four cycles on T 6 . The last 
stack is equipped with a non-trivial open string flux, to mimic the intersecting brane picture in 
type II A, as explained in section 3.4. There is also a particular mechanism of supersymmetry 
breaking through closed string 3-form fluxes. Thus, there is an explicit superpotential generated 
for the closed string moduli, which leads to an explicit dependence of the gravitino mass (1713/2), 
ps,t\u i anc j ^ e cosmological constant (V) on the moduli s, t % and u % . The cosmological constant 
is zero if the goldstino angle (8) is zero which is the same in our case. It turns out that using 
these formulas, in order for the gravitino mass to be small, the string scale is sufficiently low 
for reasonable values (0(1) in string units) of the flux. We have not assumed any particular 
mechanism of supersymmetry breaking, so we do not have an explicit expression for (1713/2) , 
ps,t\u i an( j ^y^ j n t erms Q f ^]j e moduli and have taken the string scale to be of the order of the 
unification scale. 

The model considered in [20] considers non-zero (0,3) and (3,0) form fluxes only, which leads 
to non- vanishing F s and F t% . In the T-dual version, this means that F s and F u * are non-zero, 
which is the case we examined in detail. However, in [20], an isotropic compactification is 
considered, while we allow a more general situation. 

For the calculation of soft terms, we have used the updated form of the 1/4 BPS sector 
Kahler metric as in [23], which we have also explicitly checked. In [20], the un-normalized 
general expression for calculating the soft terms has been used following whereas we use 
the normalized general expression for the soft terms in eq.([22|) [4*U] . In contrast to [20] which has 
a left-right symmetry, we have also provided an expression for the Bino mass parameter, since 
we have the SM gauge group (possibly augmented by U(l)'s) and the exact linear combination 



15 



of U(l)'s giving rise to U(1)y is known. 



5 Some phenomenological implications 

Using the formulas for the soft terms given in the previous section, we can study some aspects 
of the phenomenology of the model that has the brane setup shown in Table ^ 

Although ideally the theory should generate \i of order the soft terms and tan j3 should be 
calculated, that is not yet possible in practice as explained before. Therefore we will not specify 
the mechanism for generating the \x term. We will take tan j3 and mz as input and use EWSB 
to fix fi and b terms. 

Unlike the heterotic string models where the gauge couplings are automatically unified 1 , 
generic brane models don't have this nice feature. This is simply because in brane models 
different gauge groups originate from different branes and they do not necessarily have the 
same volume in the compactified space. Therefore to ensure gauge coupling unification at the 
scale Mx ~ 2 x 10 16 GeV, the vev of some moduli fields need to take certain values. In our 
models, the gauge couplings are determined according to eq. (|14j) . Thus the unification relations 

9 2 s =gl = \g Y ~ 0.5. (40) 

lead to three conditions on the four variables: Res and Rev, 1 where i = 1,2,3. One of the 
solutions is 

Res = 2 - 9Reui Reu 2 = 4 Rew 3 = 4. (41) 

It's interesting to note that M = 1 SUSY condition actually requires Reti2 = Rew 3 . Therefore 
although at first sight it seems that the three gauge couplings are totally unrelated in brane 
models, in this case requiring Af = 1 SUSY actually guarantees one of the gauge coupling 
unification conditions [TT] . 

After taking into account the gauge coupling unification constraint, the undetermined pa- 
rameters we are left with are m 3 / 2 , tan/5, Rett 1 , Ret 2 , Ret 3 , 0j and ji, where i = 1,2,3. The 
phases 71,2,3 enter both gaugino masses and trilinears and in general can not be rotated away, 
leading to EDMs and a variety of other interesting phenomena for colliders, dark matter, etc. 
in principle. However, for simplicity, in this paper we set all phases to be zero: 71,2,3 = 0. The 
only dependence on Ret 2 and Ret 3 is in the scalar mass squared terms and is logarithmic. For 
simplicity, we set them equal: Ret 2 = Ret 3 . Using the relation YJ |©i| 2 = 1, we can eliminate 
the magnitude of one of the 0jS but its sign is free to vary. Thus we are left with the following 
six free parameters and two signs: 

m 3 / 2 , tan/5, Rew 1 , Ret 2 , 61, 2 , sign(/i), sign(6 3 ). (42) 
1 But usually unified at a higher scale than the true GUT scale. 
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Instead of scanning the full parameter space, we show here three representative models which 
correspond to three interesting points in the parameter space. In order to reduce fine tuning, 
we require that the gluino be not too heavy. We also require that the higgs mass is about 114 
GeV, that all other experimental bounds are satisfied and that the universe is not overdosed. 
These requirements strongly constrains the free parameters. The parameters of these three 
points are shown in Table 0] 



LSP 


mz/2 


tan j3 


Rew 1 


Ret 2 


©i 


e 2 


sign(yu) 


sign(0 3 ) 


^thermal 


W 


1500 


20 


0.025 


0.01 


0.50 


0.060 


+ 




~ 


H 


2300 


30 


0.025 


0.01 


-0.75 


0.518 


+ 


+ 


~ 


B-H mixture 


2300 


30 


0.025 


0.01 


-0.75 


0.512 


+ 


+ 


- 0.23 



Table 4: Parameter choices for three particular models. All masses are in GeV. We set Ret 3 = 
Ret 2 . | ©3 1 will be fixed by the condition I] |Oj| 2 = 1. Res, Rew 2 and Rew 3 are determined by 
requiring gauge coupling unification. 

Using the values of the moduli, one can calculate the string scale Ms- It is indeed between 
the unification scale and the Planck scale. Notice that since t ~ 1/T, we are in the large radius 
limit of compactification and perturbation theory holds good. 

From the parameters shown in Table we can calculate the soft terms at high scale. 
They are shown in Table EJ We use SUSPECT |3H] to run the soft terms from the high scale 





Mi 


M 2 


M 3 


A t 


m 


W LSP 


-1288 


156 


146 


-728 


~ m 3 / 2 


H LSP 


849 


2064 


-336 


633 


~ m 3 / 2 


B-H 


866 


2040 


-336 


640 


~ m 3/2 



Table 5: Soft terms at the unification scale. The input parameters for calculating the soft terms 
are shown in TableHJ mo denotes the average of scalar masses. In both models, they are roughly 
the gravitino mass. The sign of the trilinears is according to the convention used by SUSPECT. 
It should be kept in mind that this is opposite to the convention used in supergravity formulas. 

to the weak scale and calculate the sparticle spectrum. Most string-based models that have 
been analyzed in such detail have had Bino LSPs. The three models we examine give wino, 
higgsino and mixed bino-higgsino LSP, respectively. The gluino masses 2 in the three models 
are (9(500GeV). They are significantly lighter than the gluinos in most of existing supergravity 

2 The radiative contributions to the gluino pole mass from squarks are included. 
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and superstring based MSSM models. These three models satisfy all current experimental 
constraints. If the LSP is a pure bino, usually the relic density is too large to satisfy the 
WMAP[49 constraint. For the third model, the LSP is a mixture of bino and higgsino such 
that its thermal relic density is just the measured Qcdm- For the first two models, LSPs 
annihilate quite efficiently such that their thermal relic density is negligible. Thus if only 
thermal production is taken into account, the LSP can not be the cold dark matter (CDM) 
candidate. But in general, there are non-thermal production mechanisms for the LSP, for 
example the late gravitino decay [SUj or Q-ball decays [HUES- These non-thermal production 
mechanisms have uaually been neglected but could have important effects on predicting the 
relic density of the LSP. Since at present, it is not known whether non-thermal production 
mechansims were relevant in the early universe, we have presented examples of both possibilities. 

Late gravitino decay actually can not generate enough LSPs to explain the observed flcDM 
for the first two models. This is because the gravitino decays after nucleosynthesis. Thus to 
avoid destroying the successful nucleosynthesis predictions, the gravitino should not be pro- 
duced abundantly during the reheating process. Therefore, LSPs from the gravitino decay can 
not be the dominant part of CDM. 

Another mechanism for non-thermal production is by the late Q-ball decay. Q-balls are 
non-topological solitons[HI]. At a classical level, their stability is guaranteed by some con- 
served charge, which could be associated with either a global or a local U(l) symmetry. At 
a quantum level, since the field configuration corresponding to Q-balls does not minimize the 
potential globally and in addition Q-balls by definition do not carry conserved topological num- 
bers, they will ultimately decay. Q-balls can be generated in the Affleck-Dine mechanism of 
baryogenesis|53]. Large amounts of supersymmetric scalar particles can be stored in Q-balls. 
Final Q-ball decay products will include LSPs, assuming R-parity is conserved. If Q-balls decay 
after the LSP freezes out and before nucleosynthesis, the LSP could be the CDM candidate 
and explain the observed flcDM- The relic density of N\ can be estimated as(5H| 

/0.7\ 2 / m# \ / lCT 7 GeV 2 \ /100MeV\ / 10 \ 1/2 

For our first two models, (crv) ~ lCT 7 GeV 2 . Thus if the temperature (To) when Q-balls decay 
is about 100 MeV, we will have ~ 0.22. One may attempt to relate this number to the 
baryon number of the universe since Q-balls may also carry baryon number. But probably, this 
wouldn't happen in these models at least at the perturbative level, because then baryon number 
is a global symmetry. Therefore Q-balls wouldn't carry net baryon number. They may carry 
lepton number because of lepton number violation. But since is well below the temperature 
of the electroweak (EW) phase transition, sphaleron effects are well suppressed so that the net 
lepton number cannot be transferred to baryon number. Hence, in the Q-ball scenario, baryon 
number probably is not directly related to Qt; . 
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Figure 1: Cosmic ray positron excess due to LSP annihilations. Crossed points are HEAT 
experimental data. The dotted line is from the standard cosmological model without taking 
into account the LSP annihilation contributions. The solid line is for the W LSP model, the 
dashed line is for the H LSP model and the dashed -dotted line is for the mixed B-H LSP 
model. Bs is the boost factor. 

If we assume ~ 0.23 by either thermal (for the third model) or non-thermal production 
(for the first two models), then there will be interesting experimental signatures. One of them is 
a positron excess in cosmic rays. Such an excess has been reported by the HEAT experiment [55.. 
and two other experiments (AMS and PAMELA) jSHl have been planned which will give im- 
proved results in the future. In fact the HEAT excess is currently the only reported dark matter 
signal that is consistent with other data. In both the wino LSP and the higgsino LSP models, 
LSPs annihilate to W bosons quite efficiently. In particular, W + decays to a positron(and 
a neutrino). We used DARKSUSY[o"7] to fit the HEAT experiment measurements fSHj. The 
fitted curves are shown in Figure d The "boost" factors, which characterize the local CDM 
density fluctuation against the averaged halo CDM density, in both models are not large. Small 
changes in other parameters could give boost factors very near unity. Thus both models give a 
nice explanation for the measured positron excess in the cosmic ray data. For the mixed LSP 
model, because the LSPs do not annihilate efficiently, one needs a "boost" factor of (9(100) to 
get a good fit, which may be a bit too high. 



19 



6 Conclusion 

In this paper we investigated some phenomenological implications of intersecting D-brane mod- 
els, with emphasis on dark matter. We calculated the soft SUSY breaking terms in these models 
focussing on the u moduli dominated SUSY breaking scenario in type II A, in which case the 
results do not depend on the Yukawa couplings and Wilson lines. The results depend on the 
brane wrapping numbers as well as SUSY breaking parameters. Our main result is providing 
in detail the soft-breaking Lagrangian for intersecting brane models, which provides a new set 
of soft parameters to study phenomenologically. We use a rather general parameterization of 
F-term vevs, based on [3U], which can be specialized to include the recent progress on susy 
breaking by flux compactification |2D1- Ovrr results overlap and are consistent with 

those of [20]. We applied our results to a particular intersecting brane model[T2] which gives 
an MSSM-like particle spectrum, and then selected three representative points in the parame- 
ter space with relatively light gluinos, in order to reduce fine-tuning, and calculated the weak 
scale spectrum for them. The phenomenology of the three models corresponding to the three 
points is very interesting. The LSPs have different properties. They can be either wino-like, 
higgsino-like, or a mixture of bino and higgsino. All of them can be good candidates for the 
CDM. 



7 Appendix 

Here, we list the formulas for the soft scalar mass parameters and the trilinear parameters. 
It turns out that for the matter content in Table 2, we have the following independent soft 
parameters : 

• Trilinear parameters (A): 

a QlHuUr = -v/ f m3/2 + ©s) + [(©a - © 2 ) + - *(0 - m 3 , 



47T 

+©i - - + *(£,))] siu(2vra) ] (44) 

A QL H dDR = [27r(e 2 + e 3 ) + [(e 8 - e 2 ) + ^ a 3 6 ) - *(0 - *( 



C*Cl/ 



47T 

+©i " " *(0 + *WD)1 sin (2™) ] (45) 

A Q L H d D R = A LHdER ] Aq lHuUr = A LHuNr 
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• Scalar Mass parameters (m 2 ): 



2 

m 2 QL = [-12n{Q 1 ^ (6l b ) - - (0 2 2 - 3 2 )(M>(0 2 6 ) + sm(2na) + 

3 {-20^02 - e 3 )(*'(^) - <f(0i)) + 0i 2 (^'(^) + tf^i)) + 
(0 2 - 3 ) 2 (*'(0 2 6 ) + sm 2 (27ra) + tt [(-4tt) + 

3{e?(*(&) - *(0 + (0 2 - e 8 ) a (*(^) - - 

20i(©2 - © 3 )(*(0 + *(&))} sin(4vra)] ] (46) 

2 

^ = [-12* {Gx^^) - *(0 - (©2 2 - 8s 2 )(*(0 + sin(27ra) + 

3 {-20^02 - 3 )(^'(C) " *'(0 + ©i 2 (*'(C) + *'(0 + 
(02 - © 3 ) 2 WL) + *'(0) sin 2 (2vr«) + vr [(-4vr) + 

s {e?(*(£j - + (©2 - e 8 ) 2 (*(0 - - 

20i(02 - 3 )(*(O + *(0) M*™) ] ] (47) 

2 

ml R = - n0D) - (02 2 - e 3 2 )(*(C) + *(C))} ^(27T«) + 

3 {-20 x (02 - e 3 )(*'(0 - + ©i 2 (*'(C) + %)) + 

(02 - 03) 2 (*'(#L) + sin 2 (2vra) + tt [(-4tt) + 

s {oimoU - *(0 + (© 2 - e 3 ) 2 (*(0 - *CO - 

20i(©2 - e 3 )(^(0 + sin(47ra)] ] (48) 

The above formulas are subject to the constraint |©i| 2 + | ©2 1 2 + | ©3 1 2 = 1, as can be seen below 
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